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ABSTRACT 

Characterization of the composition, surface properties, and atmospheric conditions of exoplanets is 
a rapidly progressing field as the data to study such aspects become more accessible. Bright targets, 
such as the multi-planet 55 Cancri system, allow an opportunity to achieve high signal-to-noise for the 
detection of photometric phase variations to constrain the planetary albedos. The recent discovery 
that that inner-most planet, 55 Cancri e, transits the host star introduces new prospects for studying 
this system. Here we calculate photometric phase curves at optical wavelengths for the system with 
varying assumptions for the surface and atmospheric properties of 55 Cancri e. We show that the 
large differences in geometric albedo allows one to distinguish between various surface models, that 
the scattering phase function cannot be constrained with foreseeable data, and that planet b will 
contribute significantly to the phase variation depending upon the surface of planet e. We discuss 
detection limits and how these models may be used with future instrumentation to further characterize 
these planets and distinguish between various assumptions regarding surface conditions. 
Subject headings: planetary systems - techniques: photometric - stars: individual (55 Cancri) 



1. INTRODUCTION 

The discovery of transiting exoplanets was a break- 
through in the exoplanet field. This is not simply be- 
cause it provided an additional avenue through which to 
secure their detection, but because of the vast amount 
of planetary characterization information they yield ac- 
cess to. The investment in improving photometric tech- 
niques and precision for transit surveys has led to the pe- 
ripheral consequence of detecting more subtle planetary 
signatures. Photometric phase curves of exoplanets is 
one such signature, though technically difficult to achieve 
due to the relatively low signal amplitudes. Examples 
of observed phase variations in the Infra- Red (IR) from 
Spitzer observations include HD 189733b (iKnutson et al.l 
[2009a) and HD 149026b (|Knutson et al.ll2009b[) . Exam- 
ples i n the optical includ e Kepler observations of HAT- 
P-7b (jWelsh et al.ll20ldl ) and p hase variations dete cted 
in the light curve of CoRoT-lb (|Snellen et al.ll2009f ). 

The most valuable transiting planets are the ones 
which orbit bright host stars since they lend themselves 
towards greater signal-to-noise measurements. Follow- 
up of known radial velocity (RV) pla nets is a key way 
to achieve this goal ()Kane etall[2009h . A recent exam- 
ple is that of the planets orbiting 55 Cancri (HD 75732, 
HIP 43587, HR 3522 hereafter 55 Cnc ) a br ight {V = 
5.95) G8 dwarf star. IValenti fc Fisched ([2005!) predict a 
relatively old age for th e 55 C nc system of 9.5 ± 4.4 Gyr, 
although iFischer et al.l (|2008l ) estimate a more modest 
age range of 2-8 Gyr. Thus this is an interesting system 
since the host is of similar spectral type and age to our 
own sun and yet the planetary configuration is substan- 
tially different from our own system. 

The first planet in the system was discovered by 
iButler et alJ (I1997D and the second and third were de- 
tected bv lMarcv et al.l (|2002l ). The fourth Neptune-mass 
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planet, 55 Cnc e, was found by 'McArth ur et al.l ([20041), 
with an originally deduced orbital period of 2.82 days . 
The fifth planet was discovered by iFischer et al.l ()2008l) , 
leaving us with the currently known five-planet sys- 
tem. Subsequent study o f the aliases in the power 
spectrum of the RV data by iDawson fc Fabrvckvl mW) 
found that the true period of the e planet most likely 
0.74 days ra t her th an 2.82 days. This was confirmed by 
|W inn et all ()2011l ) who successfully detected the tran- 
sit of 55 Cnc e using high-precision photometry from the 
Micr ovariabilit y and Oscillations of STars (MOST) satel- 
lite (Walk er et al. 2003). Further confirmation appeared 
through dete ction of the transit i n the IR via Spitzer 
observations (|Demorv et al.|[2"011al ). 

This system presents an opportunity to search for 
phase variations for a known system around a very bright 
star where we know much more about the relative plane- 
tary inclinations in the system than other similar multi- 
planet systen is. Phase variations for gas giants have been 
descri bed by iKane fc Gelinol (|2010[ ) and iKane fc Gelinol 
(|2011| ) and the albedos and heat r edistribution propertie s 
for such planets are discussed bv lCowan fc Agol ( 201 It ). 
Here we apply these techniques to model the phase vari- 
ations of th e 55 Cnc sy s tem u sing the revised orbital 
parameters. I Winn et"aI1 ()2011[ ) report the detection of 
phase modulations which are inconsistent with their ex- 
pectations, though more phase coverage to improve the 
amplitude measurement and account for the presence of 
the outer planets would assist in their interpretation. 
The flux ratio of a planet with radius Rp to the host 
star measured at wavelength A is defined as 



e(a. A) 



A) 
A(A) 



(1) 



where a is the phase angle of the planet. This flux ratio 
consists of three major components; the geometric albedo 
Ag{X), the phase function g{a, A), and the inverse-square 
relation to the star-planet separation r. In Section [3] we 
describe the system characteristics which are used as in- 
put into the flux ratio model for each of the planets. 
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In Section |4] we calculate the total system phase vari- 
ations based upon three different models of the inner 
planet. Finally, in Section [5l we assess the detectabil- 
ity of the individual and combined planetary signals and 
future prospects for discriminating between the surface 
models for the e planet. 

2. SCIENCE MOTIVATION 

Detection of planetary phase variations are currently 
a challenge to detect due to their relatively low ampli- 
tude. It is reasonable therefore to assess the science re- 
turn from such difficult observations. For planets which 
undergo secondary eclipses, the albedo may be deter- 
mined thro ugh careful modeling of the eclipse data (see 
for example iDeming et al.l (|2011[ )). This has the distinc- 
tive advantage of being able to measure the key planetary 
parameters of radius and orbital inclination. However, 
without complete phase variation analysis, this albedo 
determination is generally only valid at or near the sub- 
stellar point of the planet. Albedos are sensitive to such 
aspects as cloud formation depth, reflective condensates 
in the upper atmosphere, and the scattering properties 
of the surface/atmosphere. Thus there is a strong degen- 
eracy between the albedo and back-scattering properties 
when the planet is only observed at zero phase angle. By 
determining the phase function through precision obser- 
vations, one may gain further insight into these aspects 
which are inaccessible via secondary eclipse observations. 

Furthermore, continued observations of the phase func- 
tion contain information on the combined reflective prop- 
erties for all planets in the system. For planets whose 
presence and/or orbital properties are unknown, this 
can lead to ambiguity in the interpretation of the data 
(|Kane fc G elino 2010). However, in cases such as the 55 
Cnc system, the planets presence and orbital parameters 
are well determined and thus the degeneracy can be re- 
moved leading to constraints on the scattering properties 
or the outer planets. As we shall demonstrate, planet b 
contributes significantly to the total phase variation in 
this case. Even if the planets are not known to transit, 
the phase cur ves can be used to con strain the inclination 
of the orbits (|Kane fc Gelindl201lD . 

3. SYSTEM CHARACTERISTICS 

Here we describe the system characteristics which are 
used as input for the phase model, both measured and 
derived. These are summarized in TablefTJ along with the 
predicted phase variation amplitudes which are described 
in more detail in Section |4l 

3.1. Orbital Considerations 

We adopt t he complete system orbital solution shown 
in Table 10 of iDawson fc Fabrvckvl (|2010[ ). This is a self- 
consistent model which converges on the solution which 
includes the correct period for planet e of 0.74 days. The 
fit also has a smaller rms scatter of the residuals than the 
equivalent fit which reaches the old e period of 2.82 days. 
These orbital properties are shown in Table [1] including 
the period P, the mimimum mass MpSini, the semi- 
major axis a, the eccentricity e, and the argument of 
periastron uj. Note that the phase models presented here 
also account for the eccentricities present in the orbits. 
For the host star, we adopt a stellar mass of Mi, = 0.94 ± 



0.05 Mq (iFischer et al.ll2008D and a stellar radius of 
0.943 ± 0.010 Rr^ (|von Braun et alll2011| l. 



3.2. Planetary Radii 

An important property for considering the amplitude 
of the planet-to-star flux ratio from a given planet is 
the planetary radius. This quantity is normally only 
available for a planet whose transits reveal it to us, but 
can also be derived from the estimated planetary mass 
and stellar properties for the planet in qu estion. For 
planet e, we adopt the radius measured by I Winn et al.l 
(|201lD of 2.00 or 0.179 Rj. We select this ra dius 
rather than that measured by iDemorv et al.l ()2011al ) be- 
cause we are considering phase variation effects at similar 
passbands to MOST rather than Spitzer. For the more 
massive outer planets in the system , we c alculat e radii 
estima tes based upon the models of lBode nheinicr~ et al.l 
(|2003l) which take into account both the planetary mass 
and the stellar flux received at their respective semi- 
major axis. The results of these calculations are shown 
in Table □ 

3.3. Orbital Inclinations 

The orbital inclinations of the planets are unknown ex- 
cept for that of plan et e. We adopt the value of i = 90° 
for this planet from IWinn et al.l (|2011| ) . The evidence 
thus far is that none of the outer planets transit the host 
star. Here we assume that this is indeed the case and cal- 
culate the maximum inclination which satisifi es this cri- 
teria u sing the methods described in Kane fc von BraunI 
()2008f ). The inclination in this case is then given by 



cos« 



Rn + R* 



where r is the star-planet separation, given by 

a(l - e2) 



1 + e cos / 



(2) 



(3) 



and is evaluated at w + / = 7r/2. 

The above assumption for the maximum inclination 
hinges somewhat on the system orbit s being close to 
copla nar. Astrometry performed by iMcArthur et al.l 
(|2004f ) indicate that planet d may be misaligned with 
the edge-on orbit of planet e, although this is based upon 
preliminary work and the outer-most planet contributes 
a negligible amount of flux to the combined phase curve. 
The results from the Kepler mission have revealed many 
multiple-transiting system which implies that th ose sys- 
tems are remarkably coplanar (Lath am et al.l [2011). 
An excelle nt example is the six -planet system orbiting 
Kepler- 11 (ILissauer et al.ll201lD. Theoretical modeling 
performed bv iTremaiiie fc Pond (|2011D ind ependentlv 
supports the claim bv ILissauer et al. that near- 

zero mutual inclinations are favored for multi-planet sys- 
tems. Stability analyses of RV systems have also been 
used to determine coplanarity, such as the cases of the 
HD 10180 system (iLovis et al.l l201lD and GJ 876 sys- 
tem (|Bean fc Seifahrtl I2009D . " It should be noted how- 
ever that dynamical stability over long timescales can be 
achieved through interaction of giant planets or the in- 
fluence of an external perturber (jGuillochon et al.ll2011l : 
iMalmberg et al.ll2002D . 
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TABLE 1 

Planetary Orbital Parameters and Derived Characteristics 



Planet 


Pt 

(days) 


Mp sin i t 
(Mj) 


at 
(AU) 


et 


a;t 
(deg) 


iJpt 

{Rj) 


(deg) 




Flux Ratio (lO"'^) 
Rocky Molten Atmosphere 


e 


0.736537 


0.026 


0.016 


0.17 


181 


0.18 


90.0 


0.15 


3.37 20.23 


4.89 


b 


14.6507 


0.825 


0.115 


0.010 


139 


1.10 


87.5 


0.16 




3.26 


c 


44.364 


0.171 


0.240 


0.005 


252 


1.01 


88.8 


0.17 




0.70 


f 


259.8 


0.155 


0.781 


0.30 


180 


0.96 


89.6 


0.22 




0.11 


d 


5169.0 


3.82 


5.74 


0.014 


186 


1.07 


89.9 


0.50 




0.004 



t Prom Table 10 of|D awson fc Fabrvckvl 1(201^ ). 

t With the exception of planet e, based upon Bod enheimer et all 1 12003 ) models. 
* Mean geometric albedo assuming a thick atmosphere. 



3.4. Geometric Albedos 

Planetary albedos can span a very large range of values 
depending upon both the surface conditions a nd location 
of the planet. The theoretical models of Su darskv et alj 
(|2005f l show that there is a dependence of gas giant ge- 
ometric albedos on star-planet separation due to the 
removal of reflective condensates from the upper atmo- 
sphere. This was quantified by iKane fc Gelinol ()2010l ) 
who also generalized this dependence to eccentric or- 
bits. Note that this does not take into account a vari- 
able surface albedo or the thermal response of the sur- 
face/atmosphere to changing incident fiux. We use these 
models to estimate the geometric albedos for the four 
outer-most planets. 

However, the super-Earth planet e is a special case 
due to its smaller size and proximity to the host star. 
In Section 0] we consider three possible surface models 
which entail different albedos and scattering properties. 
The models are that of a rocky surface, a molten suface, 
and a thick atmosphere. The atmosphere model uses the 
geometric albedos described above with the same star- 
planet separation dependence. Shown in Table [T] are the 
mean calculated albedos for all of the planets based upon 
the thick atmosphere assumption. The rocky and molten 
surface models use the measured albedos of Mercury and 
lo as templates respectively, where the values were ex- 
tracted from the JPL HORIZONS System^. Note that 
the composition and tidal force s on lo lead to a highly 
variable surface albedo (Simonell i et al.ir2001[ ). However, 
we are considering the integrated fiux from the planet 
and so the mean albedo of lo is a useful approximation. 
Mercury's geometric albedo is 0.106 and has a density of 
5.427 gcm""^. For lo, the geometric albedo is 0.6 and the 
density is 3.530 gcm~'^. The density for 55 Cnc e is cal- 
culated from the properties in Table [T] to be 5.6 gcm^"^. 
This is comparable to that of Mercury and lo although 
the reflective properties at the surface are independent of 
the composition. Additionally, the similar bulk density 
to Mercury implies a greater concentration of volatiles 
resulting in high densities at the core with relatively low- 
density near the surface. 

3.5. Variable Incident Stellar Flux 

The geometric albedo can have a spatial dependence 
which is due to a variable amount of stellar fiux being 
received at various points on the planetary surface. Con- 
sider the cases of Mercury and Venus. Mercury experi- 
ences a large range of day-side surface temperatures due 

^ |http: //ssd. jpl .nasa.gov/?horizons| 



to a steep temperature gradient between the equator and 
the poles. The 3:2 spin-orbit resonance of Mercury re- 
sults in the periodic cooling off of the surface as it crosses 
the terminator into the night-side of the planet. In the 
case of Venus, the slow retrograde rotation would also al- 
low for extreme temperature gradients both on the day- 
side and at the day/night boundary if it were not for 
the thick atmosphere which is exceptionally efficient at 
redistributing the trapped thermal radiation. 

From the planetary radius and semi-major axis for 
55 Cnc e shown in Table [TJ the stellar ffux at the or- 
bital distance of the planetary poles is ~ 10% less than 
the flux at the equator. Spreading this flux onto the sur- 
face of the planet at a given latitude then yields a surface 
flux which is further reduced by the cosine of the angle 
the host star is from zenith at that location. However, 
we assume that the planet is tidally locked such that no 
cooling of the day-side surface ever occurs. For the var- 
ious surface models described in Section |4] we then con- 
sider a constant albedo for the day-side surface. For the 
molten surface model in particular, the variable stellar 
flux is sufficient at all latitude to result in the necessary 
sustained temperatures for melting the surface silicate 
materials. This is described further in Section |4l 

4. PHASE VARIATIONS 

Here we simulate the phase variations of the 55 Cnc 
system. The phase algori t hms a re based upon the for- 
malism of iKane fc Gelinol (|2010l ). For each planet, the 
phase angle a is defined to be zero at superior conjunc- 
tion and is described by 

cos a — sin(aj + /) sin i (4) 

where / is the true anomaly and i is the inclination of the 
orbit. The phase function of a Lambert sphere assumes 
isotropic scattering of incident flux over 27r steradians 
and is described by 

sin a + {-K — a) cos a 



(5) 



and is used for the rocky and molten surface models de- 
scribed below. For the atmosphere model, we ado pt the 
empirically derived phase function of lHiltorJ (jl992l ) which 
is based upon observations of Jupiter and Venus. This 
approach contains a correction to the planetary visual 
magnitude of the form 

Am(a) = 0.09(a/100°) + 2.39(a/100°)2 - 0.65(a/100°)^ 
leading to a phase function given by 



(7) 




Fig. 1. — The flux ratios (solid line) of the five 55 Cnc system planets to the host star, where the surface of the e planet is assumed to be 
rocky with a geometric albedo of 0.1. The time scale in each plot has been set to show one complete orbital phase of planets e (top-left), b 
(top-right), c (bottom-left), and f (bot tom -right). In each plot, the dashed line indicates the normalized phase function for the planet on 
which the figure is phased (see Section l4.1|l . These plots show that the combined phase signatures of all five planets are dominated by the 
two inner-most planets (top-left), beyond which little information can be extracted for the others without exquisite photometric precision. 



which we refer to as the Hilton phase function and ahows 
for non-isotropic (cloud) scattering. Here we confine our 
study to optical wavelengths centered on 550 nm. This 
places the study near the peak response of the Kepler and 
MOST detectors. It is possible for there to be a thermal 
component to the phase variation at thes e wavelengths, 
such a s that predicted by the models of Dcmorv ct al. 
()2011h) for Kepler- 7b. In addition, Welsh et al,, (2010) 
speculate that the Kepler observations of HAT-P-7b may 
include a thermal component due to the planet not emit- 
ting as a blackbody as is often assumed. Depending upon 
the specific surface scattering properties, this may lead 
to a non-negligible underestimate of the phase amplitude 
for the e planet that will be interesting to resolve when 
more precision photometry is acquired. Note that stud- 
ies of flux rati o dependencies on wa velength have been 
undertaken bv lSudarskv et al] (|2005f ). 

4.1. Rocky Surface Model 

Given the size, mass, and density of the planet e, as 
well as its extreme proximity to the host star, it is likely 
that any atmosphere that may once have existed has 
since evaporated and been stripped away by the stel- 
lar irradiation. We thus first consider the case that the 
surface of the planet is of a rocky form, similar to the 
surface of Mercury. The density of the planet implies a 



different composition to Mercury (indeed to all the so- 
lar terrestrial planets), but the heavier materials likely 
reside closer to the center of the planet than the surface. 

Shown in Figure [T] are the expected phase signatures 
of the 55 Cnc system using the rocky surface model for 
the planet e, where the flux ratio refers to the combined 
flux of all the planets to that of the host star. Hence the 
phase variation effects for all planets are included in each 
panel, but panels are zoomed-in on the orbital phase of 
the e, b, c, and f planets respectively and are thus phased 
on those particular planets in each case. In all four pan- 
els, the dashed line corresponds to the phase function 
of the planet on which the flgure is phased. The phase 
function as shown has been normalized to the y-scale of 
each plot to show the time-dependent contribution of the 
planet to the total phase curve. For the purposes of this 
simulation, we first assume that all planets are located at 
periastron passage, then we move time forward to start 
the phase curve where the outer-most planet in each plot 
is located at a phase angle of zero. This can be seen in 
the phase function shown in each panel which is at max- 
imum value at an orbital phase of zero. Therefore these 
simulations represent a specific orbital configuration, the 
effect of which we discuss further in Section 

The phase curves in Figure [T] show that the total phase 
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Fig. 2. — As f or the top two panels of Figurc[T] but the e planet is now assumed to have a molten surface with a geometric albedo of 0.6 
(see Section rot . 



variation of the system is dominated by planets e and b 
whose independent flux ratio amplitudes are almost iden- 
tical; 3.4 X 10-6 and 3.3 x 10"^ respectively. The flux 
ratio amplitude of the c planet is almost an order of mag- 
nitude less; 7.0 x IQ-^. Note the asymmetric modulation 
for the e and f planets (top-left and bottom-right panels) 
due to their respective eccentric orbits. The maximum 
flux ratio amplitude of the f planet is 1.1 x 10"^ and 
occurs near the periastron passage near phase 0.7. The 
outer-most planet, d, is sufficiently far away from the 
host star that it contributes a negligible amount of flux 
to the total flux ratio; 4.0 x 10"^. These results are sum- 
marized in Table [T] 

4.2. Molten Surface Model 

|W mn "eTall (I20T1I) estimate a surface temperature of 
2800 K at the substellar point if the planet is tidally- 
locked. Even if the heat is somehow redistributed, the 
equilibrium temperature will still be as high as 1980 K. 
This is well above the melting po int for igneous rock. 
In add ition, the best-fit solution bv lDawson fc FabrvckyI 
()2010[ ) contains a non-zero eccentricity for planet e 
which would create a "super-Io" effect, as described by 
I Barnes et al.l (|2010l l. We thus here consider the entirely 
plausible case of a molten surface for the inner planet. 

Figure [2] is equivalent to the top two panels of Figure 
[1] in that they are zoomed-in to the phases of planets e 
and b. The much higher albedo causes the e planet to 
become dominant in the phase curve shown in the right 
panel (note the different ordinate scales). The e planet 
now has a flux ratio amplitude of 2.0 x 10"^ which is 
almost an order of magnitude higher than that for the b 
planet. 

4.3. Atmosphere Model 

The e planet is unlikely to harbor an atmos phere under 
the e xtreme conditions of its environment (jWinn et al.l 
l2011t ) . Here we consider this possibility for completeness 
and as a direct comparison with the other two presented 
scenarios for the surface of planet e. In this case we use 
the Hilton rather than the Lambert phase function to 
represent the non-isotropic scattering of the atmosphere. 
The non-uniform incident stellar flux described in Sec- 
tion 13.51 is insufficient to allow an non-uniform albedo 



since the reflective consendates will be removed from the 
upper atmosphere regardless of latitude under such ex- 
treme temperatures. As shown in Table [l] the mean ge- 
ometric albedo for the atmosphere model is only slightly 
higher than that for the rocky surface due to these re- 
flective condensates being effectively removed. 

In Figure [3] we see the combined calculated flux ratio 
zoomed to the phases of the e and b planets where, once 
again, the ordinate scales have increased relative to Fig- 
ure [1] This model produces a flux ratio amplitude for 
planet e of 4.9 x 10"^ which slightly exceeds the ampli- 
tude expected from the b planet. We discuss being able 
to distinguish between this model and the rocky surface 
model in Section [S) 

4.4. Orbital Configurations and Phase Modulation 

The simulation results provided above describe a spe- 
cific starting configuration for the planets which is based 
upon their periastron passages and phase angles. How- 
ever, for observations at some random epoch, this config- 
uration will be arbitrary in nature and the phase curves 
shown here will not necessarily match that which is ob- 
served. In particular, the relative phases of the two plan- 
ets which contribute the bulk of the total reflected light, 
e and b, will result in phase modulations. Thus there 
will be optimal orbital configurations for which to search 
for the phase signatures. The maximum amplitude of 
the flux ratio will approximately occur where both of 
the planets are simultaneously close to zero phase an- 
gle. However, the extraction of only the phase informa- 
tion for planet e will be best achieved when planet b is 
near inferior conjunction. Fortunately, the orbital peri- 
ods of these two planets are relatively short which allows 
the scheduling of such optimal observating times to be 
straightforward . 

If one requires both high signal-to-noise throughout the 
measurement, and folding of multiple orbits of planet e, 
then one needs to be mindful of the contributions of the 
outer planets as they progress in their orbits. For each 
orbit of planet e, the outer planets will shift in phase by 
0.05, 0.017, 0.0028, and 0.00014 respectively. Through 
a single orbit this will have minor effects on the total 
phase variation, but over multiple orbits, will begin to 
show modulation effects which will dampen the signa- 
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Fig. 3. — As for the top two panels of Figure[T] but the e planet is now assumed to have an atmosphere with non-isotropic scattering 
properties (see Section 14.31 1. 



ture of planet e, thus resuhing in an incorrect estimate of 
the flux ratio amplitude. Deriving accurate ephemerides 
through further RV measurements will help to curtail 
such effects in photometric monitoring of the phase vari- 
ations. 

5. SIGNAL DETECTABILITY 

Here we assess the detectability of the planetary flux 
ratios and the possible effects of stellar variability. 

5.1. Instrumentation Requirements 

From the expected phase amplitude of the two inner- 
most planets of 55 Cnc, the instrumentation requirement 
for successful detection of the phase signatures is pho- 
tometry with a relative photometric precision of at least 
r^ 10^^. For the e and b planetary phase variations, 
long-term stability of high-precision photometry is not 
required. The Kepler mission is already achieving this 
precision for significantly fainter stars, though it should 
be noted that the exquisite photometers of Kepler were 
designed to perform such a task. An example of the high- 
standard of Kepler precis ion is the det e ction of phase and 
ellipsoidal variations by I Welsh et al.l ()201Clt ) using only 
the Ql Kepler data, where the amplitude of this variation 
is3.7xl0"^. Using the STIS ins trument on the Hubble 
Space Telescope (HST). [Brown et al. (2001) obtained a 
precision of 1.1 x 10"** per 60 second integration observa- 
tion during primary transits of HD 209458b. Such high 
time resolution is not required for phase variation ob- 
servations and so binning these data would improve the 
rms scatter. A future mission that would allow such ob- 
servations of 55 Cnc to be carried out include the James 
Webb Space Telescope (JWST) with the NIRCam instru- 
ment, though wavelength range of this instrument would 
include a substantial thermal component of the phase 
variation which would need to be accounted for. The 
minimum precision mentioned above would adequately 
confirm or rule out a molten surface for planet e. Further 
RV characterization of the orbits for the 55 Cnc planets 
will allow one to accurately predict both the amplitude 
of the predicted phase signature and times of maximum 
and minimum flux ratios. This knowledge will help to 
distinguish the phase signatures from instrumental drift 
effects. 



Future possabilities also exist from the ground, al- 
though one needs to also contend with the offsets from 
night-to-night variations. There are several large tele- 
scopes under development that are capable of meeting 
the challange of very high photometric precision, such as 
the European Extremely Large Telescope (E-ELT), the 
Thirty Meter Telescope (TMT), and the Giant Magel- 
lan Telescope (GM T). It has also been demonstrated by 
iColon et alJ (|20lol ) that precision photometry of < 0.05% 
can be achieved with the 10.4m Gran Telescopio Canarias 
through the use of narrow-band filters. These was con- 
ducted for observing the signatures of known transiting 
planets which is possible to achieve within a single night. 
Longer term monitoring of phase signatures will require 
careful accounting for the aforementioned nightly varia- 
tions in addition to the airmass corrections throughout 
a night. 

5.2. Stellar Variability of 55 Cnc 

If one is able to accomplish the required level of pho- 
tometric precision, the greatest impediment to studying 
the planetary phase variations will be the i ntrinsic stel- 
lar var iability. An analysis of Kepler data bv lCiardi et al.l 
(|2011h found that most dwarf stars are stable down to 
the the precision of the Kepler spacecraft, with G-dwarfs 
being the most stable of the studied spectral types. The 
main cause of photometric variability in F-G-K stars is 
starspots and rotation. The rotation period for 55 Cnc 
has been measured on numerous occasion s through pho- 
tometric variations. ISimpson et al.' (201(f) calculate a ro- 
tation period of 44.1 days and Fischer et al. (2008) mea- 
sure a rotation period of 44.7 days. iWinn et al.l (j201ll ) 
also observed variation of the order 10~^ which is as- 
sumed to be due to both stellar activity and rotation. 
For the c planet where the orbital period is close to the 
rotation period of the star, the variation due to phase and 
rotation may be difficult to disentangle. The peaks in the 
power spectrum from a fourier analysis of the photometry 
may separate to a degree where the starspot variability 
can be isolated from the phase signature. The known 
phase of the planet from the RV analysis will be the 
greatest aid in discriminating these two signals. It should 
also be noted that there is an M dwarf binary companion, 
55 Cnc B, with an angular separation of 84.7"(1150 AU) 
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(|Mugrauer et al.1 [20061 ) . so it is unlikely to be inside a 
photometric aperture. 

6. CONCLUSIONS 

The multiple-planet system of 55 Cnc presents many 
opportunities to understand and characterize a system 
around a solar analog with very different characteristics 
to our own system. The discovery that the inner-planet 
transits enhances the opportunities since we can place 
greater constraints on experiments designs to investigate 
these properties. Here we have specifically addressed the 
method of detecting the phase variations of the planets. 
These results show that the inner two planets have flux 
ratio amplitudes which arc comparable to what has al- 
ready been detected by Kepler around much fainter stars 
and will be accessible to next-generation ground and 
space-based observing platforms if not before. Which of 
the two planets dominates the phase signature depends 
upon if the planet has a rocky or molten surface, with the 
molten surface model producing a dominance of planet 
e to the signature due to the higher reflective proper- 
ties of the surface. The phase variation of the e planet 
for a rocky surface is almost indistinguishable from that 
of an atmosphere model at the level of the precision re- 
quirements for detection, 10^^. These two models have 
very similar surface albedos and it makes little difference 
whether one assumes isotropic or non-isotropic (cloud) 
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scattering at such small star-planet separations. The 
outer planets do not present a significant impediment 
for the detection of the e and b phase variations but one 
needs to be aware of the expected phase modulations 
if multiple phase observations are undertaken to boost 
signal-to-noise. We have not considered here the thermal 
component of the flux from the e planet, whose variation 
could be sign ificant if ind e ed th e planet is tidally-lock as 
suggested bv lWinn et all (|2011h . 

Even though the transit depth is relatively small, the 
high signal-to-noise possible from transiting planets such 
as 55 Cnc e demonstrate the value of such objects. The 
best hope then for more such transiting planets around 
bright stars lies through investigation of known exoplan - 
ets discovered using the RV technique ()Kane "erall[2009l) . 
Projects such as the Transit Ephemeris Refinement and 
Monitoring Survey (TERMS) seek to accomplish just 
that which will hopefully lead to further characterization 
opportunities. 
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